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We present results ol QCD corrections to Higgs boson production at the CERN Large Hadron Col- 
lider. Potentially large logarithmic contributions from multiple soft -gluon emission are resummed 
up to all order in the strong coupling. Various kinematical distributions, including the Higgs trans- 
verse momentum, are predicted at the O(a^) level. Comparison is made to outputs of the popular 
Monte Carlo event generator PYTHIA. 
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I. INTRODUCTION 

One of the fundamental questions of the Standard 
Model (SM) of elementary particle physics is the dynam- 
ics of the electroweak symmetry breaking. Within the 
SM, the Higgs mechanism postulates the existence of a 
scalar field, the elementary excitation of which is called 
the Higgs boson. Four experimental collaborations at 
the LEP II collider search for the Higgs boson in the 
e + e~ — > Z°H process up to 202 GeV center of mass en- 
ergy. DELPHI and L3 set a preliminary exclusion limit of 
itlh > 96 GeV on the Higgs mass, followed closely by the 
limit of m H > 91 GeV set by ALEPH and OPAL §. Ac- 
cording to recent preliminary information, the combined 
lower limit is close to m# ~ 106 GeV 0. Global fits 
to electroweak observables appear to prefer a low mass 
Higgs particle, with a mean value close to 100 GeV, and 
less than 250 GeV within 95 percent of confidence HQ . 

Among other aims the main goal of the next proton- 
proton accelerator, the 14 TeV center of mass energy 
Large Hadron Collider (LHC) at CERN, is to establish 
the existence of the Higgs boson and to measure its ba- 
sic properties. At the LHC a light SM Higgs boson will 
be mainly produced through the partonic subprocess gg 
(via top quark loop) — > H X || . It can be detected, after 
about 1.5 years of running with a statistical significance 
of at least 4, in its H — > 77 decay mode, if its mass is in 
the 100-150 GeV range Q. If the Higgs mass is higher 
than about 150 GeV, then its H -> Z°W Z° decay mode 
is the cleanest and most significant 0. In this letter 
our focus is on the Higgs boson production, and in our 
numerical illustration we choose ran = 150 GeV. 

According to earlier studies, a statistical significance 
on the order of 5-10 can be reached for the inclusive 
H — > 77 and for the H +jct — > 77 + jet signals, although 
actual values depend on luminosity and background esti- 
mates. In Ref. it was found that in order to optimize 



the significance it is necessary to impose a 30 GeV cut on 
the transverse momentum of the jet, or equivalently (at 
next-to-leading order precision), on the transverse mo- 
mentum (Qt) of the photon pair. With this cut in place, 
extraction of the signal in the Higgs + jet mode requires 
the precise knowledge of both the signal and background 
distributions in the medium to high Qt region. 

To reliably predict the Qt distribution of Higgs bosons 
at the LHC, especially for the low to medium Qt region 
where the bulk of the rate is, the effects of the multi- 
ple soft-gluon emission have to be included. One ap- 
proach to achieve this is parton showering || . Although 
the universality of this method makes it a very power- 
ful tool, present drawbacks of this ansatz are the lack of 
the proper normalization (which takes into account the 
full fixed order QCD corrections), the lack of exact ma- 
trix elements even in the high Qt region, and the lack 
of uniqueness of the prediction ("tunability" ). There is 
ongoing work to correct these problems . 

A more reliable prediction of the Higgs Qt can be 
obtained utilizing the Collins-Soper-Sterman (CSS) re- 
summation formalism |Ic|Jll| , which takes into account 
the effects of the multiple soft-gluon emission while re- 
producing the rate, systematically including the higher 
order corrections. It is possible to smoothly match the 
CSS result to the fixed order one in the medium to high 
Qt region, thus obtaining the best prediction in the full 
Qt region fl2|| . Compared to fixed orders, the resummed 
result depends on a few extra parameters. These pa- 
rameters are new renormalization scales (C*j, only two of 
which are independent) ]Tl| , and a few universal, non- 
perturbative parameters (gi), which are extracted from 
present experiments and then used to predict the results 
of future ones p3[ . In this letter, we use this formalism 
to calculate the total cross sections and Qt distributions 
of Higgs bosons at the LHC. 

Our results here, together with the resummed calcu- 
lations for the diphoton and Z° boson pair production 



-|17|, provide a consistent set of QCD calculations of 
the transverse momentum (and other) distribution(s) of 
the Higgs bosons and their backgrounds at the LHC. 
These results systematically include both the multiple 
soft-gluon effects and the finite order QCD corrections, 
and can be used to tune the shower MC's which experi- 
mentalists extensively use when extracting the Higgs sig- 
nal, or can be utilized independently by means of the 
RcsBos Monte Carlo event generator [fl2|. 



II. ANALYTICAL RESULTS 

Within the SM ||,|l8| , as well as in the minimal super- 
symmetric standard model (MSSM) with small tan (3 [0 , 
the dominant production mode of neutral Higgs bosons at 
the LHC is gluon fusion via a heavy quark loop. The low- 
est order cross section of this process is formally 0(a 2 ) 
in the strong coupling. Fixed order QCD corrections to 
this production mechanism are known to substantially 
increase the rate. For a light Higgs boson the O(a^) to 
0(a 2 ) if -factor is in the order of 2 (cf. Fig. Q). The full 
O(a^) calculation is not completed yet, but the real emis- 
sion |ll| and the virtual contributions |2(J are separately 
available. Resummed calculations, taking into account 
the soft-gluon effect, were also performed to estimate 
the size of the uncalculated higher order corrections [^l) , 
as well as to predict the shape of the Higgs transverse 



momentum distribution |22 23y . 

In this work we use the Collins-Soper-Sterman (CSS) 
soft-gluon resummation formalism to calculate the QCD 
corrections from the multiple-soft gluon emission. Cal- 
culations similar to this were earlier performed in Refs. 
p2| , p3| . Our present calculation improves these by in- 
cluding 0(a^) terms in the Sudakov exponent, by ap- 
plying the state of the art matching to the latest fixed 
order distributions, by using a QCD improved gluon- 
Higgs effective coupling |g3| , by utilizing an improved 
non-pcrturbativc function, and by including the effect of 
the Higgs width. 

We also utilize the approximation that the object 
which couples the gluons to the Higgs (the top quark 
in the SM) is much heavier than the Higgs itself. This 
approximation is not essential to our calculation and can 
be released by including the complete Wilson coefficients 
with all the relevant masses. The heavy quark approxi- 
mation in the SM has been shown to be reliable within 5 
percent for ran < 2m t ]2q- |27|] , and still reasonable even 
in the range of m# ~ 2m t |2l|]. It has also been shown 
that the approximation remains valid for the Qt distri- 
bution in the large Qt region, provided that m# < m% 
and Qt < nit pl|. In this work we assume that the 
approximation is valid in the whole Qt region. In the 
MSSM the heavy quark approximation is also a reliable 
ansatz for the case of the light Higgs boson and small 
tan/3 when the Yukawa coupling of the bottom quark 



is negligible (cf. |2l[] and references therein). Using the 
CSS formalism we resum large logs of the type \u(Q/Qt) 
in the low Qt region, and we match the resummed result 
to the fixed order calculation which is valid for high Qt 
]l2[ . We also include the qg and qq subprocesses which, 
in combination, can constitute up to 10 percent of the 
total rate, depending on the Higgs mass [^5) . 

The resummed differential cross section of a neutral 
Higgs boson, denoted by </>° in the SM or MSSM, pro- 
duced in hadronic collisions is written as 
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The kinematical variables Q, Qt, and y are the invariant 
mass, transverse momentum, and rapidity of the Higgs 
boson, respectively, in the laboratory frame. The parton 
momentum fractions are defined as x\ — e y 'Mt / 'v~S \ and 
x 2 = e- v M T /\/S, with M T = ^Q 2 + Q T , and y/S being 
the center-of-mass (CM) energy of the hadrons hi and 
1%2. The lowest order cross section is 
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where Gf is the Fermi constant, and k^, the QCD cor- 
rected effective coupling of the Higgs boson to gluons in 
the heavy top quark limit (cf. Ref. [pTJ ) , is defined as 
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where a s is the strong coupling constant in the MS 
scheme with 5 active flavors, and mt denotes the pole 
mass of the top quark. 

The renormalization group invariant kernel of the 
Fourier integral W gg (b^, Q, x\, X2, 61,2,3), an d the Qt 
regular term Y{Qt, Q, x\,X2, C4), together with the vari- 
ables 6* and C± to C4 , are given in Ref. |^3| . The defini- 
tion of W gg , contains the Sudakov exponent 



S(Q,K,C ll C 2 ) = 
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In the perturbative expansion of the A (a s (/l), C\) and 
B (a s (fl), Ci, C2) functions we follow the notation of Ref. 
p9[ . In our present calculation, we include the process 
independent next-to-next-to-leading order coefficient 
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in the expansion of the A(a s {ji),C\) function, where 
C a — 3 is the Casimir of the adjoint representation of 
SU(3), No = 3 is the number of colors, and Nf = 5 is 
the number of active quark flavors. With the inclusion of 
the only missing next-to-next-to-leading order con- 
tribution in the Sudakov exponent is the term, which 

is suppressed by 1/ln (gr) with respect to A^ 2 \ and by 

a s with respect to B^\ This is illustrated by the expan- 
sion of the asymptotic part of the cross section: 
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where i and j label incoming partons. While the A^ 
coefficient contributes to the above series via 
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the B^ coefficient only occurs as 
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where our notation coincides with that of Ref. p9[ . 
Hence, the contribution from the uncalculated B^ 2 \ com- 
pared to that from A^ , is expected to be smaller because 
of the additional log weighting the A^> coefficient. To 
estimate the size of the contribution from B^ 2 \ we fol- 
low the usual practice in a perturbative calculation by 
varying the renormalization constants (C± and C2) in the 
Sudakov factor by a factor of 2. The results are shown 
in Fig. |. _ 

The form of our non-perturbative function W^ g p coin- 
cides with the one used for the gg — > 77 process in Ref. 
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where the Casimir of the fundamental SU(3) represen- 
tation is Cp = 4/3. The values of the non-perturbative 
parameters gi are defined in Ref. (To]. The uncertainties 



of the resummed distribution, stemming from the non- 
perturbative function, were found to be in the order of 
5 percent (cf. poj ). In the high Qt region Eq. (|l|) is 
matched to the fixed order perturbative result (at the 
0(a s )) of Ref. in the manner described in Ref. [pj. 



III. NUMERICAL RESULTS 

The analytic results are coded in the ResBos Monte 
Carlo event generator jl2 17 1, which uses the following 
electroweak input parameters Wm : 



G F = 1.16639 x 10 -t) GeV , m z = 91.187 GeV, 
m w = 80.36 GeV. 



(10) 



As in the background calculation |T|] , we use the canon- 
ical choice of the renormalization constants (C\ = C3 = 
2e~~< E = C and C 2 = C 4 = 1 0), the NLO expressions 
for the running electromagnetic and strong couplings 
a(n) and ag(jj,), as well as the NLO parton distribution 
function set CTEQ4M (defined in the modified minimal 
subtraction, MS, scheme) |p2]| . We set the renormaliza- 
tion scale equal to the factorization scale: fiR — /if = Q- 
In the choice of the non-perturbative parameters, we fol- 
low Ref. [H| . Since we are not concerned with the decays 
of Higgs bosons in this work, we do not impose any kine- 
matic cuts. We defer the more extensive study, including 
various decay modes and QCD backgrounds, to a future 
publication Ref. [Q. 

Fig. |l| displays Higgs boson production cross sections 
via the gluon fusion process at the LHC, calculated with 
various QCD corrections in the SM as the function of the 
Higgs mass. The ratio of the fixed order O(o\) (dashed) 
and the lowest order 0(a 2 ) (dotted) curves varies be- 
tween 2.0 and 2.3. We note that less than 2 percent of the 
0(a^) corrections come from the qg and qq initial states 
for Higgs masses below 200 GeV. The resummed curve is 
about 10 percent higher than the 0(a^) one, as expected 
based on the findings that the CSS formalism preserves 
the fixed order rate within the error of the matching (the 
latter being higher order) |L2|j , The resummed rate is 
close to the 0(a\), because we used the 0(p? s ) fixed or- 
der results to derive the Wilson coefficients which are 
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a resummed cal- 



utilized in our calculation. In Ref. 
culation estimated the size of the 0(a s ) corrections, and 
a typical value of 1.5 of the O(af) to 0(af) if-factor can 
be inferred from that work. Based on this, we also plot 
the O(a^) curve rescaled by 1.5, to illustrate the possible 
size of the O(af) corrections and to establish the nor- 
malization of our resummed calculation among the fixed 
order results. 

Fig. compares the Higgs boson transverse momen- 
tum distributions calculated by ResBos (curves) and by 
PYTHIA H (histograms from version 6.122). The mid- 
dle solid curve is calculated using the canonical choice 
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FIG. 1. SM Higgs boson production cross sections at 
the LHC via gluon fusion as the function of the Higgs mass, 
with QCD corrections calculated by soft -gluon resummation 
(solid), at fixed order O(al) (dashed), and without QCD cor- 
rections at O(o.g) (dotted). The 0{o? a ) curve is scaled by 1.5 
(dash-dotted, c.f. Ref. ||2l| ) to estimate the 0{a i B ) result. 



for the renormalization constants in the Sudakov expo- 
nent: C\ — Co, and C% = 1. To estimate the size of 
the uncalculated term, we varied these renormaliza- 
tion constants multiplying both by 1/2 and 2. The upper 
solid curve shows the result for Cl = Co/2, C2 = 1/2, and 
the lower solid curve for Cl = 2Co,C2 = 2. The band 
between these two curves gives the order of the uncer- 
tainty following from the exclusion of B^ 2 \ The typical 
size of this uncertainty, e.g. around the peak region, is in 
the order of ±10 percent. The corresponding uncertainty 
in the total cross section is also in the same order. The 
dashed PYTHIA histogram is plotted without altering its 
output. The normalization of PYTHIA, as that of any 
parton shower MCs, is the lowest order 0(a 2 ) for this 
process. The default PYTHIA histogram is also plotted 
after the rate is multiplied by the factor K = 2 (dotted) . 
The shape of the PYTHIA histogram agrees reasonably 
with the resummed curve in the low and intermediate Qt 
(< 125 GeV) region. For large Q T the PYTHIA predic- 
tion falls under the ResBos curve, since ResBos mostly 
uses the exact fixed order 0(a^) matrix elements in that 
region (c.f. Ref. jl^l), while PYTHIA still relies on the 
multi-parton radiation ansatz. PYTHIA can be tuned to 
agree with ResBos in the high Qt region (dash-dotted), 
by changing the maximal virtuality a parton can acquire 
in the course of the shower, i.e. the Q^ al parameter, 
from the default value to the partonic center of mass en- 
ergy s. In that case, however, the low Qt region will 
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FIG. 2. Higgs boson transverse momentum distributions 
calculated by ResBos (curves) and PYTHIA (histograms). 
The default (middle) ResBos curve was calculated with the 
canonical choice of the renormalization constants, and the 
other two with doubled (lower curve) and halved (upper 
curve) values of Ci and Ci. For PYTHIA we show the origi- 
nal output with default input parameters (dashed) , the same 
rescaled by a factor of K = 2 (dotted), and a curve calculated 
by the altered input parameter value Q max = s (dash-dotted). 
The lower portion, with a logarithmic scale, also shows the 
high Qt region. 



have disagreement, since the normalization in PYTHIA 
is conserved, so events in the low Qt region are depleted. 



IV. CONCLUSIONS 

In this letter we presented Higgs boson production 
rates and Qt distributions for the LHC, including O(al) 
fixed order QCD and multiple soft-gluonic corrections by 
means of the CSS resummation formalism. We showed 
that the resummed rate recovers the fixed order 0(a^) 
rate, as expected within the CSS formalism. We investi- 
gated the uncertainty of the resummed prediction due to 
uncalculated terms in the Sudakov exponent. We found 
that the shape of the resummed prediction in the low Qt 
region is in reasonable agreement with the default result 
of PYTHIA. 
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